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PERTURBATIONS IN THE CURVATON MODEL 
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Recently, the primordial non-Gaussianity in the curvaton model has been predicted assuming 
sudden decay of the curvaton. We extend the calculation to non-instantaneous decay by 
employing SN formalism. The difference between the sudden-decay approximation and our 
numerical result is larger than 1% only if the non-linearity parameter is small, —1.16 < /nl < 
60. Thus it is safe to use the sudden-decay approximation when deriving constraints for the 
curvaton model from WMAP3 (/nl < 114), but with the Planck forecast |/nl| < 5 one should 
employ the fully numerical result. Often, the curvaton perturbations Sa have been assumed 
O |- to be small compared to the background value of the curvaton field ao. Consequently, the 

variance A 2 = (Sa 2 ) /a 2 , has been assumed to be negligible. However, the measurements of 
CMB or large-scale structure perturbation amplitude do not constrain the variance if the main 
■ contribution to it comes from the ultraviolet (UV) scales, i.e., from smaller than observable 

' scales. We discuss how, even in this case, observational constraints on non-Gaussianity set an 

upper bound to the small scale variance, A^y < 90. 
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Most inflationary models give rise to nearly Gaussian primordial curvature perturbation. Typ- 
ically, prediction for the non-linearity parameter /nl in single-field models is of the order of e 
- the slow roll parameter which must be < 10 _1 to guarantee the near scale-invariance of the 
primordial perturbations. In principle, measurement of /nl would give valuable information on 
Q-f the inflaton potential, but unfortunately such a tiny non-Gaussianity is likely to remain unob- 

servable. The current upper bound from the WMAP three- year data^ is |/nl| < 114 while 
£j . Planck is expected to bring this down to |/nl| ^5 5, which is still orders of magnitude larger 

than the typical inflationary prediction. 

Nevertheless, there are classes of multi-field models that can lead to an observable non- 
Gaussianity. One well-motivated example is the curvaton mo deP In addition to the inflaton 6 
there would be another, weakly coupled, light scalar field (e.g., MSSM flat direction), curvaton 
(7, which was completely subdominant during inflation so that the inflaton drove the expansion 
of the universe. The potential could be as simple as^ V = ^M 2 <p 2 + ^m 2 a 2 . At Hubble exit 
both fields acquire some classical perturbations that freeze in. However, the observed cosmic 
microwave (CMB) and large-scale structure (LSS) perturbations can result from the curvaton 
instead of the inflaton, if the inflaton perturbations are much smaller than 10~ 5 . To simplify the 
analysis, in this talk a , we assume that the curvature perturbation from inflaton is completely 
negligible G, <C Cm an d the curvature perturbation from curvaton C a is such that it leads to the 
observed amplitude of perturbations. 

After the end of inflation the inflaton decays into ultra-relativistic particles ( "radiation" ) the 
curvaton energy density still being subdominant. At this stage the curvaton carries pure entropy 
(isocurvature) perturbation instead of the usual adiabatic perturbation. Namely, the entropy 
perturbation between radiation and curvaton is S ra = 3( C — C a ) « — 3( a . Since the observations 
have ruled out pure isocurvature primordial perturbation- , a mechanism — curvaton decay 
into "radiation" — that converts the isocurvature perturbation to the adiabatic one is needed 
at some stage of the evolution before the primordial nucleosynthesis. 
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As the Hubble rate, H, decreases with time, eventually H 2 < m 2 , and the curvaton starts to 
oscillate about the minimum of its potential. Then it behaves like pressureless dust ( "matter" , 
p a oc a~ 3 ) so that its relative energy density starts to grow with respect to radiation (p r oc a -4 ). 
Finally, the curvaton decays into ultra-relativistic particles leading to the standard radiation 
dominated adiabatic primordial perturbations. 6 However, this mechanism may create from the 
initially Gaussian curvaton field perturbation a strongly non-Gaussian primordial curvature 
perturbation. The more subdominant the curvaton is during its decay the more non-Gaussianity 
results in. Since the time of the decay depends on the model parameters (such as the curvaton 
mass m and decay rate T a ), the observational upper bounds on non-Gaussianity provide a 
method to constrain these parameters. 

In the simplest case the (possible) non-Gaussianity results from the second-order correction 
to the linear result 

C = Cl + k 2 = Cl + |/NLCl 2 - (i) 
2 5 

Here £i is proportional to the Gaussian field perturbation at the Hubble exit, so it is Gaussian, 
and Ci is X 2 distributed. In this talk we briefly derive /nl in the sudden-decay approximation 
using a slightly different and more general approach than fn^SI anc [ then, for the first time, 
present the results in the case of non-instantaneous decay of curvaton. For the full derivation 
and discussion of our results se where we, in addition to the second-order calculation, go to 
the third order, and finally derive fully non-linear results and full probability density function 
(pdf) of the primordial C, in the long- wavelength limit. Since, in the early universe, all today's 
observable scales are super-Hubble, we take advantage of the separate universe assumption 
throughout the calculations. 



2 Sudden-decay approximation 

In the absence of interactions, fluids with a barotropic equation of state, such as radiation 
(Pr = Pr/3) or the non-relativistic curvaton (P a = 0), have a conserved curvature perturbation 

( i (t,af) = 6N(t,x) + ~ P * (2) 

where 5N = N(t, x) — N(t) with N(t, x) being the perturbed (i.e. local) number of e-folds 
of expansion until time t and N(t) being the average expansion. At the first order this fully 
non-linear definition reduces to the usual definition: Cn = —tfti — with ib\ = —b\N. 

__ Pi 

Applying (J2J) for the curvaton during its oscillation but before the decay, we have 



1 

Ca{t,x) = - In 



Pa(t,x) 
Pa(t) 



(3) 

SN=0 



where p a (t,x) = \m 2 a 2 {x) is evaluated on spatially flat (5N = 0) hypersurface, and o~(x) it 
the local amplitude of oscillation. The time of the beginning of the curvaton oscillation, ij n , 
is defined in terms of the local Hubble rate as H(ti n ,x) = m. Since H 2 = (87rG/3)p to t, the 
constant time t = i; n surface is a uniform-total density hypersurface. 

In general, <7i n depends non-linearly on the field value at Hubble exit a*. Thus we writ!3 
0"in = g{o~*) = g(cf*) +g'5a* + ^g"(5a*) 2 + . . ., where ' = d/da*. (For exactly quadratic potential 



^Hadn't we assumed negligible inflaton curvature perturbation, Q, « 0, some "residual" isocurvature would 
have resulted if the curvaton was sub-dominant during its decay. This would have lead to an interesting mixture 
of correlated adiabatic and isocurvature perturbations which was studied in^. Following the guidelines of^ 
our calculation should be straightforward to generalise. It should be noted that observations do not rule out a 
correlated isocurvature component if it is less than 20% of the total primordial perturbation amplitude!^ 
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g" and higher derivatives vani sh, b ut even a slight deviation from quadratic potential can change 
the resulting /nl considerably'^ via g" '.) Substituting this into © and expanding up to second 

order we obtain Cr = Cri + K ct2 + • • • witn Cri = f^ "* and Ct2 = — § (l - gg"/g' 2 ) Clii where 
<5<7* is well described by a Gaussian random field a s we a ssume the inflaton and curvaton to be 
uncoupled (or only weakly coupled), see e.g. 

RefdPEEl We can express this result in terms 
of the effective non-linearity parameter for the curvaton perturbation, analogous to Eq. (|T]l. 
/jJ L = — | ^1 — ^j/jrj- Hence we find /nl = —5/4 for the curvaton ( a in the absence of any 
non- linear evolution (g" = 0). If the curvaton comes to dominate the total energy density in the 
universe before it decays, so that C = Cn then this is the generic prediction for the primordial 
/nl in the curvaton model, as emphasised by^. 

Assume now that the curvaton decays instantaneously at time idee (before it has become 
completely dominant) on a uniform-total density hypersurface corresponding to H ~ T a , i.e., 
when the local Hubble rate equals the decay rate for the curvaton (assumed constant). Hence 

Pr(tdec,x) + Pa(tdcc,x) = p(t d ec) , (4) 

where we use a bar to denote the homogeneous, unperturbed quantity. Note that from Eq. (j2J we 
have C = oN on the decay surface, and we can interpret £ as the perturbed expansion, or U 5N" . 
Assuming all the curvaton decay products are relativistic (P = p/3), we have that £ is conserved 
after the curvaton decay. The local curvaton and radiation densities on this decay surface may be 
inhomogeneous. Indeed we have from Eq. (J2J ( r = C + \ l n (Pr/f>r) and Q a = ( + g In (p a /p a ) or, 
equivalently, p r = p r e 4 ^ T ~^ and p a = p a e 3 ^~'''. Requiring that the total density is uniform on 
the decay surface, Eq. (j3J), then gives the simple relation (1 — ^l fT ,dec) e ^ < ' r ~^ + ^decC 3 ^ 17- ^ = 1 > 
where fi CT ,dcc = p~a/(p r + At) Idee is the dimensionless density parameter for the curvaton at the 
decay time. This equation can be rewritten in the form 

e 3C = 3±r e 3C + ^3 e -C ; (5) 

where r = ^CZc ' Recallin g Ec l- © the LHS of Ec l- © is P*(tm> *)/ P*(tm) = V^*]- 
As Eq. Q is a fourth degree equation for e^, the primordial curvature perturbation, £, as 
a function of Gaussian <r* can be solved exactly. Remarkably, Eq. (J5J) was derived using fully 
non-linear definitions. Hence we have found an exact fully non-linear solution for £ (as opposite 
to only second-order calculations in the literature). Expanding this solution in 5a* up to second 
order we then find /nl, up to third order the so called ^nl, etc. Since the exact solution is quite 
long, it turns out to be easier to directly expand RHS of Eq. © up to any wanted order, and 
then equate it to the LHS, i.e., to g 2 [a*(x)]/g 2 [a 1 „] order by order. 

Up to second order the solution is £i = rfc-l = r |^"<5o"* and (2 = [-§p (l+W/*/ 2 ) — 2 — r] 
These give the non-linearity parameter (^Q) in the sudden-decay approximation^^! 

/NL = |:(l + ^)-f-f • (6) 

In the limit r — > 1, when the curvaton dominates the total energy density before it decays, we 
recover the non-linearity parameter /^ L of the curvaton /nl — * /nl = — f (l — 9d" /d' 2 ) ■ On the 
other hand we may get a large non-Gaussianity (|/nl| S> 1) in the limit r — > 0, where we have 
fNL^i{l+gg"/g' 2 )- 

3 Non-instantaneous decay 

After the curvaton has decayed the universe is dominated by "radiation", with equation of 
state P = p/3, and hence the curvature perturbation is non-linearly conserved on large scales. 
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Figure 1: (a) The non-linearity parameter /nl as a function of curvature perturbation transfer efficiency r — 
Cr,out/Co-,m- The analytical approximative, i.e., sudden-decay result (red dashed line) crosses zero at r — 0.58 and 
is negative for r > 0.58. The exact numerical result (black solid line) crosses zero at r — 0.53. (b) r vs pi n . 



Although the sudden-decay approximation gives a good intuitive derivation of both the linear 
curvature perturbation and the non-linearity parameter arising from second-order effects, it is 
only approximate since it assumes the curvaton is not interacting with the radiation, and hence 
£ CT remains constant on large scales, right up until curvaton decays. In practice the curvaton 
energy density is continually decaying once the curvaton begins oscillating until finally (when 
T > H) its density becomes negligible, and during this process Co- does evolvi^. Another 
problem with results derived from the sudden-decay approximation is that the final amplitude 
of the primordial curvature perturbation, and its non-linearity, are given in terms of the density 
of the curvaton at the decay time which is not simply related to the initial curvaton density, 
especially as the decay time, H ~ r CT , itself is somewhat ambiguous. 

In the non-instantaneous case we can still define the first order transfer efficiency "r" of 
the initial curvaton perturbation to the output radiation perturbation, but now it must be 
calculated numerically from the definition r = Cir,out/Ci<r,in- It turns out to be a function solely 
of the parameter^ pi„ = [0 CT (^) 1/2 ]in = m^i^ 2 , where M|j = 1/(8ttG). 

We start the numerical integration of the background Friedmann eqn together with radiation 
and curvaton continuity eqns (with ±T a p a as a source term) from uniform-density surface at ti n 
and end it when curvaton has completely decayed at some suitable uniform-density hypersurface 
H(t en d,x) = -ffend- Then the fully non-linear primordial curvature perturbation will be ( = 
Cr,out = 5N(t en d, x). Repeating the calculation with different initial values a m we find the 
function N{a- m ). From this we then calculate the first and second order quantities by employing 
5N = N'5o* + \N"{5a*) 2 + . . . = Ci + \C2 + ■ ■ ■, and finally 8 / NL = ^N" / (N') 2 . From Fig. 1(a) 
we see that if /nl > 60 (r < 0.02) or /nl < —1-16 (r > 0.95), the sudden-decay result differs 
from the non-instantaneous decay result less than 1%. C In the both cases, as r — > 1 we have 
/nl —* — f ■ As r — > 0, the sudden-decay result is /nl — ► f ^(1 + 99" Id' 2 ) ~ 1-67, whereas 
the non-instantaneous decay gives /nl - * \h (1 + gg" /g' 2 ) — 2.27. We can write the full result 
in a form /nl = |f(l + gg" ' / 'g' 2 ) + §/i(r)/r 2 , where the function h(r) is defined by equation 
r' = [2r + h(r)]g' /g, and determined numerically in the non-instantaneous case. For the sudden 
decay we have from © h(r) = —^r 2 — |r 3 . 



c After this talk the non-instantaneous decay calculation was done in™™ employing second order perturbation 
equations instead of SN formalism which we use. The results agree after taking into account that the authors of 

compare /nls resulting from the fixed pi n while we prefer to compare /nl when the sudden decay and non- 
instantaneous decay produce same r, or, in other words, the same observable first order curvature perturbation. 
To achieve this, we need to start from different pi n for sudden decay than for non-instantaneous decay. The 
non-trivial mapping between p ln and r is demonstrated in Fig. 1(b) for the non-instantaneous decay. 
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4 Large variance on small scales 



Finally, we consider the possibility of a large small-scale variance**-^ A 2 = {(^§ hL ) 2 )- Let us 
name the observable CMB scales as infrared (IR) and smaller scales as ultraviolet (UV) so 
that \ir S> Xjjv 3> . As the observations require ((f) ir < 10~ 9 , they set constraint for 
|r 2 A| J? . However, observations do not directly constrain A 2 ^. Allowing large A^y we find 
/nl = f(l + A 2 / y)^(l + <7</ / /s ,/2 ) + f h(r)/r 2 , where h(r) remains same as without large variance. 
Thus, large homogeneous small-scale variance modifies the first term of /nl only. Recalling that 
r < 1, and -54 < / NL < 114 (from WMAP30), we find an upper bound l + A 2 ^ < § x 114 = 91. 
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